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We report the emergence of large-scale hyperuniformity in microfluidic emulsions. Upon peri¬ 
odic driving confined emulsions undergo a first-order transition from a reversible to an irreversible 
dynamics. We evidence that this dynamical transition is accompanied by structural changes at all 
scales yielding macroscopic yet finite hyperuniform structures. Numerical simulations are performed 
to single out the very ingredients responsible for the suppression of density fluctuations. We show 
that as opposed to equilibrium systems the long-range nature of the hydrodynamic interactions 
are not required for the formation of hyperuniform patterns, thereby suggesting a robust relation 
between reversibility and hyperuniformity which should hold in a broad class of periodically driven 
materials. 

PACS numbers: 05.65.-l-b,47.57.Be, 05.70.Ln 


What is the most effective way to homogeneously fill 
space with an ensemble of particles? At thermal equi¬ 
librium, an obvious effective strategy would be to endow 
the particles with interactions promoting the formation 
of a crystal. In a d-dimensional system, thermal fluc¬ 
tuations would spontaneously organize the particles into 
an ordered state where the fluctuations of the number 
of particles enclosed in a box of size i would scale as 
ANf ~ crystals are hyperuniform [l|. At large 

scales, they are much more homogeneous than a random 
set of points with number fluctuations of the order of 
the box volume £‘^. However c^^stals are not the only 
patterns being hyperuniform [l|, . Over the last decade 

much attention has been devoted to disordered structures 
displaying miniature density fluctuations. As it turns out 
such hyperuniform patterns have been shown to display 
outstanding optical properties such as complete photon¬ 
ics band gaps Until very recently the only two 

controlled strategies to engineer hyperuniform materials 
were based on numerical optimization techniques iii, 


or the jamming of athermal hard spheres [7l-lll|. In 2015, 
two sets of numerical simulations have demonstrated 
that hyperuniformity emerges when ensembles of parti¬ 
cles driven out of equilibrium approach a critical absorb¬ 
ing phase transition ^ 131. However, these numerical 
models based on elegant toy models lack a truly analo¬ 
gous physical system, in which hyperuniformity emerges 
from genuine physical interactions. 

In this letter we demonstrate experimentally that, 
when periodically driven, microfluidic emulsions self- 
organize into macroscopic (yet finite) hyperuniform 
structures. We evidence that maximal hyperuniformity 
is reached at the onset of reversibility of the droplet dy¬ 
namics, even though the emulsion does not display any 
critical feature and does not reach a genuine absorbing 
state either. We also identify the minimal ingredients 
required to produce hyperuniform emulsions by means 



FIG. 1: Close-up of the emulsion at the beginning of the 500th 
cycle. The active droplets are shown as open orange circles, 
the passive droplets as filled blue circles, (a) For a driving am¬ 
plitude A/A* = 0.65, the dynamics is reversible, (b) For 
a driving amplitude A/A* = 1.54, the dynamics is not re¬ 
versible anymore. Note also the markedly different structure 
in the two cases: for low amplitudes the structure is more 
homogeneous. 


of numerical simulations. Surprisingly, unlike in equi¬ 
librium systems, we show that long-range interactions 
are not necessary to yield extended hyperuniform struc¬ 
tures. We therefore conjecture that hyperuniformity is 
intimately related to reversibility in periodically driven 
systems, and should therefore be achieved in a broad class 
of hard and soft condensed matter materials. 


The experimental setup is the one used in 1^ con¬ 
sisting of flowing a monodisperse emulsion (area frac¬ 
tion: 0.36) in a microfluidic channel (0.5 cmx5 cmx27 ± 
0.1 ^m). The droplets have a diameter a = 25.5 ± 0.5 ^tm 
comparable to the height of the channel and therefore un¬ 
dergo two-dimensional motion. The inlet of the channel 
is connected to a syringe pump that drives the suspension 
sinusoidally. The emulsion is prepared in a reproducible 
initial state reached after a sequence of 10 high-amplitude 
oscillations. Then, a sequence of 10^ cycles at the desired 
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FIG. 2: (a) Fraction of active particles fa in the steady state at various oscillation amplitudes. At A = A* reversibility 
abruptly breaks down. Inset: Variations of the mean-squared strobed displacement in the flow direction {5x^) plotted as a 
function of A. (b) Centered and normalized probability density distribution of the strobed displacements of the droplets, Sx, 
in the flow direction for the experiments (top curves) and numerics (bottom curves). The numerical curves have been shifted 
for sake of clarity. The color indicates the mean fraction of active particles. See colorbar in (a), (c) Structure factor S{k) at 
various amplitudes (Experiments: thin lines, simulations: Thick lines. Same color code as in (b)). For sake of clarity the curves 
corresponding to A < A* have been shifted down by a constant value. The dashed line is a guide and corresponds to 


amplitude is applied to ensure that the measurements are 
performed in a statistically stationary state. The mean 
displacement of the droplets occurs along the main flow 
direction and is sinusoidal, and its amplitude A is the 
sole control parameter of the experiments. Data is col¬ 
lected by tracking the instantaneous position of 3 x 10^ 
droplets remaining in the field of view throughout the en¬ 
tire flow cycle. Two snapshots of the emulsion are shown 
in Fig. [T] and correspond to A = 18.3a and A = 43.3a 
respectively. 

Following [i^ the macroscopic reversibility of the sys¬ 
tem is measured by determining the fraction of active 
particles fa, which is the fraction of droplets that behave 
irreversibly. A droplet is here defined to be active if it 
does not return at the end of a cycle within the spatial 
extent of the Voronoi' cell it occupied at the start of the 
cycle. In these experiments Re <C 1, the fluid flows are 
therefore reversible in time 1^. However, above a driv¬ 
ing amplitude A*/a = 28.1 ± 0.3 the droplet dynamics 
abruptly becomes irreversible. A macroscopic fraction of 
the droplets remains endlessly active upon periodic driv¬ 
ing as illustrated in Fig. 1, and quantified in Fig. [5^ 
where fa is plotted as a function of driving amplitude. 

As the order parameter fa is based on a metric-free 
criteria, it is both affected by the changes in the dynam¬ 
ics, and in the structure of the emulsion. We now dis¬ 
entangle and elucidate these two concomitant collective 
phenomena. Let us begin with the dynamical arrest 
of the strobed dynamics. The reversible states where 
droplets retrace their steps back to their initial Voronoi' 
cell, A < A*, do not correspond to interaction-free con¬ 
formations: the passive droplets continuously interact 
with their neighbors via hydrodynamic interactions in 


the course of the cycles. They also a priori experience 
a number of weak but irreversible perturbations such as 
short-range potential interactions, or minute shape defor¬ 
mations which cannot be experimentally measured. As a 
result, even in the reversible regime, they return only on 
average to their initial position after a cycle. The transi¬ 
tion to an irreversible regime where all the particles are 
active is associated to a discontinuous amplification of 
the mean-square displacement of the strobed dynamics, 
thereby causing the escape of the droplets from their ini¬ 
tial Vorono'i cell, see Figure [2^ inset. We characterize 
the fluctuations of the strobed dynamics by the distribu¬ 
tion of the displacements at the end of each cycle in Fig¬ 
ure [2)3. At the onset of irreversibility this distribution is 
not merely widened, the statistics of the droplet displace¬ 
ments becomes qualitatively different at A* thereby con¬ 
firming the discontinuous nature of the dynamical tran¬ 
sition. Surprisingly, whereas the strobed-displacement 
statistics is Gaussian for high A, it is exponential in 
the reversible regime, echoing the existence of rare large- 
amplitude jumps when the droplets exit their Vorono’i 
cell. These intermittent displacements are akin to the 
cage jumps found in glass forming liquids [l^. Impor¬ 
tantly, these results imply that the dynamical arrest of 
the strobed dynamics does not belon g to the absorbing- 
phase-transition scenario reported in 113, [HE [13. 

We now focus on the central results of this letter and 
demonstrate the emergence of hyperuniformity. Having a 
quick look back at Fig.[l] we see that the configuration in 
the reversible regime looks much more homogeneous than 
in the irreversible regime. In order to quantify this ap¬ 
parent structural change, we compute the structure fac¬ 
tor S{k) of the emulsion at various amplitudes. Fig. j^):. 
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FIG. 3: (a), (b) and (c) Normalized density fluctuations for experiments (a), simulations with long-range interactions (b), and 
short-range interactions (c). At high amplitudes (red curves) the density fluctuations are normal, whereas in the reversible 
regime, at low amplitudes (blue curves), the density fluctuations are suppressed, (d) Fitted power-law exponents A from the 
curves plotted in (a), (b) and (c). The lines are guides to the eye. Density fluctuations are only suppressed for A < A*, (e) 
Extent of the hyperuniform regions tnu- Filled symbols: experiments, open symbols: simulations. The lines are guides to the 
eye. The system is only hyperuniform up to length-scales of ~ 10a (resp. ~ 15a) in the experiments (resp. in the simulations). 


The existence of a structural transition at the onset of 
reversibility is very noticeable from the sudden change of 
S{k) across all scales. The shift of of the high-fc peak 
of the structure factor corresponds to the change of the 
liquid-like structure reported in 1^. However the most 
striking feature at the transition occurs at low k. For 
high driving amplitudes (upper curves), S(k) plateaus to 
a finite value as k goes to 0, indicating the absence of 
long-range order. In contrast, at low amplitudes (lower 
curves) S{k) decays algebraically as k goes to 0, large- 
scale density fluctuations are suppressed [l|. The emul¬ 
sion self-organizes into a hyperuniform state. 


In order to quantify the degree of hyperuniformity and 
the extent of the hyperuniform regions, we go back to real 
space and directly measure the statistics of the droplet 
number Ni in a ixi box Figure[3K shows the vari¬ 

ations of the variance ANg = (7V|) — (Ng) normalized by 
(N/) as a function of i. Any ensemble of particles with no 
spatial correlation whatsoever has a variance that scales 
as ANg (Ne). Therefore, ANf/ (Ni) is a decreasing 
function of i for hyperuniform systems. In Figure [3^ we 
quantify the level of hyperuniformity of the emulsions by 
normalizing the variance ANf by that of a random set 
of points of the same size and density. This procedure 
is used to minimize any statistical artifacts due to the 
finite sample size. For sake of clarity we have shifted the 
curves starting from small driving amplitudes at the top. 
For large drivings (A ^ A*), A7V^/(7V^) does not show 
any significant variations. The density fluctuations are 
normal. Conversely, as the dynamics becomes reversible 
(A < A*), the emulsion becomes locally hyperuniform. 
ANgjNi indeed first decays algebraically with £ up to 


a box size £hu above which it increases. We expect the 
variations to then saturate for high values of ^ that are 
not accessible with our experimental setup. This non¬ 
monotonic behavior is very similar to the one reported 
first in 12| , at the onset of an absorbing phase transition. 


More quantitatively, the system homogeneity is quanti¬ 
fied by the exponent A > 0 defined as: ANg / {N() oc 
For a random set of points A = 0, while A = 1 for perfect 
crystals. In Fig. [31(b) the exponent A fitted for small I 
is shown as a function of the driving amplitude A. In¬ 
terestingly, the closer the system to the reversible tran¬ 
sition, the more the density fluctuations are suppressed. 
At A = A*, we find that A ^ 0.5 which is again close to 
the one reported for systems close to an absorbing phase 
transition (A 0.45) in 12| and in [l^ at intermediate 
scales. 

The extent of the hyperuniform regions is characterized 
by measuring the length £hu where AN} / (Ng) deviates 
from a decreasing power law. Fig. ISKc) [l8|. This length- 
scale undergoes non-monotonic variations with the driv¬ 
ing amplitude. When increasing A, £hu decays from ten 
to five droplet diameters around A/A* ^ 0.7. Then ap¬ 
proaching the reversibility transition £hu increases again 
to its maximal value (^ 10a) before dropping down to 
zero above A*. Two comments are in order. Firstly, we 
do not see any sign of a divergence of £hu as A approaches 
A*, which means that the emulsion never self-organizes 
into a fully hyperuniform state where the density fluc¬ 
tuations would be suppressed at the entire system scale. 
However, the typical extent of the hyperuniform regions 
are much larger than the typical distance below which 
the droplets display translational order. The pair cor¬ 
relation function of the emulsion decays exponentially 
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over distances that are at most of the order of a cou¬ 


ple of particle diameters [18|. Secondly, the rather com¬ 
plex variations of £hu contrasts with that of all the other 
structural, and dynamical quantities which only display 
a significant change at the transition point A*. A po¬ 
tential explanation is that the intrinsic slowing down of 
the strobed dynamics below A* makes the hyperuniform 
self-organization too slow to be experimentally achieved, 
although all the other (local) observables have reached a 
steady state. 

What causes this emulsion to self-organize into hyper¬ 
uniform large-scale structures? The droplets in the ex¬ 
perimental system interact through various forces: hy¬ 
drodynamic forces that are time reversible, but also 
short-range irreversible forces such as depletion, van der 
Waals, and electrostatic forces which are specific to the 
nature of the fluids and surfactants forming the emul¬ 
sion. To find out which of these ingredients are rele¬ 
vant to achieve hyperuniformity, we perform numerical 
simulations using a model containing only minimal hy¬ 
drodynamic interactions and steric repulsion due to the 
finite size of the droplets. The flow induced by a mov¬ 
ing droplet in a geometry as used in the experiment is 
described by a potential-flow dipole 0 . Each of the 
896 droplets is advected by the local flow u(r) with 
the friction being modeled through a mobility coefficient 
0 < /J, < 1. The equation of motion for droplet i is 
Ti = ^u(ri). Assuming pairwise additive interactions, 
the flow at the location of particle i is u(ri) is the sum of 
the contributions from the driving flow Uq, and from the 
flow induced uj by all other particles i 7 ^ j in the system. 
The full equations of motion are: 


ri = fj, \ Uo{t)+ 


3^i 


27r|ri. 


( 1 ) 


where (Jj is the dipole-vector associated with the particle 
j. The strength of the dipole is proportional to the ve¬ 
locity of this particle relative to the ambient fluid. The 
boundary conditions in the flow-direction are periodic, 
whereas the flow is bounded by walls in the transverse di¬ 
rection. The infinite number of dipole-images that arise 
due to this are modeled analogous as in | 20 ll and is de¬ 
scribed in the supplementary materials [l 8 |. 

This minimal model correctly accounts both for the 
dynamical and structural transitions. A reversible-to- 
irreversible transition occurs as the driving amplitude A 
is varied. Fig. [2^, open symbols. The dipole strength 
which is our sole free parameter, is set to match the ex¬ 
perimental value of A*. Without any additional adjust¬ 
ment, we see that the same dynamical changes occur as 
in the experiments. Fig. [ 2 ] 3 . The shape of the probabil¬ 
ity density abruptly changes from a Gaussian to an ex¬ 
ponential across the reversibility transition as observed 
in the experiments. Similarly, the structure above and 
below the transition are markedly different at all scales. 


Fig. [ 2 j 3 . Again the computed structure factor share the 
same salient features as the experimental one. Impor¬ 
tantly, our simulations quantitatively captures the emer¬ 
gent hyperuniform structures in the reversible regime. 
Fig. [31 Both the exponent A and the extent Ihu of 
the hyperuniform regions are consistent with the experi¬ 
mental measurements. Fig. | 35 ,c. This agreement unam¬ 
biguously demonstrates (i) that hyperuniformity chiefly 
stems from the combination of reversible hydrodynamic 
interactions and short-range repulsion, and (ii) that this 
phenomenology is robust to the very details of the inter¬ 
actions between the droplets and of their near-held flows. 

In order to gain more physical insight, we now ques¬ 
tion the importance of the long-range r~^ in 2D) na¬ 
ture of the interactions, which are very specihc to hy¬ 
drodynamics. Long-ranged interactions can yield hy¬ 
peruniformity in systems at thermal equilibrium such 
as one-component plasmas [2l|. A natural question is 
therefore: does the emergence of hyperuniformity depend 
on the long-range nature of the particle-particle interac¬ 
tions in this non-equilibrium system as well? To answer 
this question, we perform the same simulations as above, 
but apply a very short-ranged exponential screening to 
the hydrodynamic interactions of the form exp(— |rij|/a), 
keeping all the other parameters unchanged. Choosing 
a screening length of one particle diameter prevents the 
droplets from interacting over the observed hyperunifor¬ 
mity length scale while preserving the reversible nature of 
the microscopic dynamics. As is evident from Fig.|2K, the 
same reversible-to-irreversible transition, yet smoother, 
is observed thereby further demonstrating the robustness 
of our main findings. Counterintuitively, the extent of the 
hyperuniform regions is clearly not reduced by screening, 
conversely it is extended up to the entire simulation win¬ 
dow. As opposed to equilibrium systems, long-ranged 
interactions impair hyperuniform self-organization. This 
counterintuitive observation might be explained by a 
critical reversibility transition when the interactions are 
short-range, as suggested by the sharp increase of the 
structural relaxation time at A* [l^. This picture is 
consistent with the results reported in [l5j where hydro- 
dynamic interactions are short-ranged and the transition 
critical. However, a thorough finite-size scaling analysis 
would be required to unambiguously confirm that long- 
range interactions suppress criticality in our system. 

Together with that of Bill, our experimental and 
numerical results strongly suggest that any ensemble of 
particles at the onset of a reversible-to-irreversible tran¬ 
sition self-organizes into hyperuniform patterns. A broad 
class of materials are therefore expected to self-assemble 
into hyperuniform structures upon periodic driving, from 


to 


colloidal suspensions 15, 22-24|, to soft glasses [2 
shaken grains to vortices in superconductors [27 1 . 
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